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Abstract 

We report on a three-year program of coordinated X-ray and optical monitoring of 
the narrow- line Seyfert 1 galaxy NGC 4051. The rapid continuum variations ob- 
served in the X-ray spectra are not detected in the optical, although the X-ray and 
optical continuum fluxes are correlated on time scales of many weeks and longer. 
Variations in the flux of the broad H/3 line are found to lag behind the optical con- 
tinuum variations by approximately 6 days (with an uncertainty of 2-3 days), and 
combining this with the line width yields a virial mass estimate of ~ 1.1 x 10 6 Mq, at 
the very low end of the distribution of AGN masses measured by line reverberation. 
Strong variability of HeiiA4686 is also detected, and the response time measured 
is similar to that of H/3, but with a much larger uncertainty. The He n A4686 line is 
almost five times broader than H/3, and it is strongly blueward asymmetric, as are 
the high-ionization UV lines recorded in archive spectra of NGC 4051. The data are 
consistent with the Balmer lines arising in a low-inclination (nearly face-on) disk- 
like configuration, and the high-ionization lines arising in an outflowing wind, of 
which we observe preferentially the near side. During the third year of monitoring, 
both the X-ray continuum and the He n A4686 line went into extremely low states, 
although the optical continuum and the H/3 broad line were both still present and 
variable. We suggest that the inner part of the accretion disk may have gone into 
an advection-dominated state, yielding little radiation from the hotter inner disk. 
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1 Introduction 



Beginning in early 1996, the International AGN WatchF] undertook a pro- 
gram of contemporaneous X-ray and optical spectroscopic monitoring of the 
galaxy NGC 4051, the only NLS1 galaxy in Seyfert's ( fLlf) original list of high 
surface-brightness galaxies with strong emission lines. The X-ray variability 
characteristics of NGC 4051 are typical of the NLS1 class (|5|; §). X-ray ob- 
servations were made with the Rossi X-Ray Timing Explorer (RXTE), and 
optical spectra were obtained with the Perkins 1.8-m telescope at Lowell Ob- 
servatory and the 1.6-m Tillinghast Reflector on Mt. Hopkins. The purpose of 
this program has been twofold: 

(1) To determine the nature of the relationship between the X-ray and UV- 
optical continuum variations. This is a particularly interesting question 
in the case of NGC 4051 since the X-ray flux dropped to an extremely 
low level towards the end of this campaign (|17D . 

(2) To determine the size of the broad-line region (BLR) and virial mass of 
the central object via reverberation techniques (|IJ [12]). 

The emission-line cross-correlation lag r can be taken to be the light-travel 
time across the BLR, so the BLR size is given by r = cr. By combining this 
with the emission-line width V FWUM , the mass of the central source can be 
inferred to be 

M — fV ™™ CT , (1) 

where / is a factor of order unity that depends on the still unknown geometry 
and kinematics of the BLR. There is an implicit assumption that the gravita- 
tional force of the central object dominates the kinematics of the BLR; this is 
formally unproven, but at least in the case of the well-studied Seyfert 1 galaxy 
NGC 5548, the reverberation-mapping data are consistent with the required 
Vfwhm oc r -1 / 2 relationship ([14]). Virial mass estimates based on reverberation- 
mapping data are now available for nearly 40 AGNs (|T8| ; ^). 



2 Results 

In Fig. 1, we show the continuum and emission-line light curves obtained 
during this three-year program. While there is a clear lack of correlated short 

1 All International AGN Watch papers and data are available at the AGN Watch 
website at URL fittp: / /www. astronomy.ohio-state.edu/ ~agnwatch / , 
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time-scale behavior of the X-ray and optical continua, the light curves in 
Fig. 1 suggest that a correlation on longer time scales is possible. To test 
this quantitatively, we have suppressed the rapid variations, smoothing both 
the optical continuum and X-ray light curves with a rectangular function of 
width 30 days, similar to what was done in a comparison of X-ray and optical 
variability in the Seyfert 1 galaxy NGC 3516 (f7|). Cross-correlation of the 
overlapping parts of these light curves yields a lag of the optical variations 
relative to the X-ray of r = 6^^ 2 days with a correlation coefficient r max = 
0.74, i.e., the mean X-ray and optical fluxes are indeed correlated once the 
high-frequency variability is suppressed. The lag between variations in the 
two wavebands is highly uncertain, but consistent with zero or any small time 
lag expected in the continuum emitting region. 



The time delay between continuum and emission-line variations has been de- 
termined by cross-correlation of the light curves. The period from JD2450183 
to JD2450262 is particularly well-suited to this type of analysis because the 
light curves are well-sampled and the character of the variations permits an 
accurate cross-correlation measurement. The time delays, or lags, we measure 
are r = 5.9t2!c) days for H/3 and r = 4. 51^0 days for HellA4686, where the 
uncertainties are based on a model-independent Monte Carlo method (|13D. 



By combining this lag with the Doppler width of the emission line, we can es- 
timate a virial mass, as in eq. (1); for consistency with previous work (|18"| ; |]). 
we take / = 3/ a/2. Since the broad emission-line features are comprised of 
a number of different components (or contaminants), it is desirable to mea- 
sure the Doppler width of only the variable part of the emission line. In order 
to isolate the variable part of the emission line and exclude constant com- 
ponents (such as contamination from the narrow-line region), we measure 
the relevant line widths in the root-mean-square (rms) spectrum constructed 
from the spectra used to determine the time lag; the principal virtue of the 
rms spectrum is that it isolates the features in the spectrum that are actu- 
ally varying. We measure widths V FWHM = 1110 ± 190 km s _1 for B./3, and 
Vfwhm = 5430 ± 510 km s _1 for HellA4686; the latter is also strongly blue- 
ward asymmetric. Examination of archival IUE spectra show that the UV 
high-ionization lines are similarly asymmetric. On the basis of the H/? vari- 
ations, a virial mass of l.ljlo'l x 10 6 M & is inferred; unfortunately however, 
the extremely large uncertainty in the He II A4686 lag renders the virial mass 
obtained from it not very enlightening (i.e., M = 19. x 1O 6 M ), but it 
is consistent with the FL3 result. It is also important to keep in mind that be- 
cause of the unknown geometry and kinematics of the BLR, the virial mass is 
reliable to only about an order of magnitude, i.e., the systematic uncertainties 
are much larger than the errors quoted here. 
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Fig. 1. Continuum and emission-line light curves of NGC 4051 between January 
1996 and June 1998. Top panel is 2-10 keV count rate. Second panel is optical 
continuum flux in units of 10 -15 ergs s _1 cm~ 2 A -1 . Third and fourth panels are 
HP and HeiiA4686 fluxes, respectively, in units of 10 -13 ergs s _1 cm -2 . 



3 Discussion 



3.1 The Continuum 



It has already been pointed out based on these same RXTE data (|17D that the 
X-ray continuum of NGC 4051 virtually "turned off" in early 1998. However, 
the optical spectroscopic data show that the optical continuum and emission 
lines (and therefore, by inference, the ionizing UV continuum) did not disap- 
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Fig. 2. Top panel shows the X-ray light curve as in Fig. 1. Bottom left panel 
shows the optical rms spectrum based on data obtained between JD2450183 and 
JD2450262. Bottom right panel shows the optical rms spectrum based on data ob- 
tained between JD2450810 and JD2451022. 



pear at the same time. Figure 2a shows the RXTE light curve. The two lower 
panels of Fig. 2 show the optical rms spectra based on data obtained during 
the intervals marked on Fig. 2a. Figure 2b shows the rms spectrum obtained 
during a period of relatively active X-ray and optical continuum variability. 
The rms spectrum in Fig. 2c is based on spectra obtained during the period 
when the X-rays were about an order of magnitude fainter than normal. Figure 
2b shows that when the X-ray flux was high and variable, the optical contin- 
uum and the H/3 and HellA4686 emission lines were also variable. Figure 2c, 
however, shows that when the X-rays were in a faint state, both the optical 
continuum and the H/3 line were still varying, but the He II A4686 emission line 
was no longer varying. A possible interpretation of the behavior of NGC 4051 
is that the inner X-ray producing part of the accretion disk has entered an 
advection-dominated accretion-flow (ADAF) state, in which radiation is emit- 
ted with very low efficiency ( |T0t H). The outer part of the disk, which produces 
the UV-optical continuum and drives the emission lines, remains relatively un- 
affected, however. The implication is that there is a transition radius inside of 
which the disk is an ADAF and outside of which it radiates efficiently, perhaps 
like a classical thin disk (|l~6D, and the persistence of the optical continuum and 
the H/3 emission line suggests that this transition radius is somewhere between 
the regions that are most responsible for the soft X-rays and the H-ionizing 
continuum. 
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3.2 The Virial Mass and Implications for NLSls 

Reverberation-based size estimates for the broad emission lines and resulting 
virial mass estimates provide a potential means of distinguishing among the 
various NLS1 models. In Fig. 3a, we show the relationship between the BLR 
radius as measured from the Balmer-line lags as a function of the optical 
continuum luminosity for all AGNs with Balmer-line lags known to reasonable 
accuracy (JT^; ^). This compilation contains six additional AGNs that could 
be classified as NLSls as they meet the criterion Vf WH m ~ 2000 km s _1 . The 
best-fit regression line (-Rblr oc 2y°- 62±0 02 ), based on all objects (excluding 
NGC 4051), is shown as a dotted line. NGC 4051 lies approximately 2.8a 
above this regression line, although all the other narrow-line objects clearly 
fall in the locus defined by the AGNs with broader lines. Given the large 
dispersion in this relationship, reflected in the high xt( = 15.7) of the fit, the 
offset of NGC 4051 is not unusual. 

Figure 3b shows the mass-luminosity relationship for these AGNs. We show 
(a) the best-fit regression line based on all objects except the seven narrow-line 
objects and (b) that based on the narrow-line objects alone. These two fits 
are separated by typically an order of magnitude in black-hole mass; the black 
holes in the narrow-line objects are about a factor of 10 lower than those of 
other AGNs of comparable luminosity. 

How well do these results allow us to distinguish among the various explana- 
tions for the NLS1 phenomenon? We consider the possibilities: 

(1) Do the BLRs of NLSls have anomalously large radii? The position of 
NGC 4051 in Fig. 3a might suggest that this is possible, but the distri- 
bution of other narrow-line objects does not support this. Furthermore, 
as noted above, the scatter in the BLR-radius luminosity relationship is 
very large, and NGC 4051 is in a statistical sense not the largest outlier in 
this relationship (simply because other sources have smaller uncertainties 
in their measured lags). 

(2) Are NLSls simply low-inclination systems? If the BLR is a flattened 
system, at low inclination (i.e., nearly face-on) the line widths will be 
decreased by a factor sinz, but the measured emission-line lags will be 
relatively unaffected. On the other hand, assuming that the UV-optical 
continuum arises in an accretion disk at the same inclination, the appar- 
ent UV-optical luminosity is higher at lower inclination. Thus, relative to 
similar sources at intermediate inclination, the masses of low-inclination 
sources will be underestimated, and their luminosities will be overesti- 
mated, displacing the narrow-line objects in Fig. 3b towards the lower 
right. This is generally consistent with the location of all seven of the 
narrow-line objects, including NGC 4051. The line transfer function for 
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Fig. 3. (a). The relationship between size of the Balmer-line emitting region and 
optical luminosity. The dotted line is a least-squares fit to all the data, excluding 
NGC 4051. (b). The relationship between virial mass and optical luminosity. The 
dashed line is a least-squares fit to the narrow-line objects only, and the dotted 
line is a least-squares fit to all of the other points. In both panels, the filled circles 
are Seyfert 1 galaxies and the open circles are QSOs. The large triangles are those 
AGNs in which H/3 has a half width less than 2000 km s _1 . NGC 4051 is shown as 
a filled diamond. 



H/9 would provide a more definitive test of this hypothesis since it would 
allow determination of the inclination of the system. This would require 
more and better data than we have obtained in this experiment. 
(3) Are NLSls undermassive systems with relatively high accretion rates? 
Again, the distribution of the narrow-line objects, including NGC 4051, 
in Fig. 3b is consistent with this hypothesis. The narrow-line sources on 
this plot lie below the mass-luminosity relationship for other AGNs, at 
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the lower end of the envelope around this relationship. 



In summary, the hypothesis that NLSls have unusually distant BLRs for their 
luminosity is probably not viable in general, although it could apply to the 
specific case of NGC 4051. At the present time, however, we cannot distinguish 
between the low-inclination and low-mass, high accretion-rate hypotheses on 
the basis of the reverberation results alone, although considerations based on 
X-ray data favor the latter. Indeed, it is entirely possible that both effects 
(i.e., low inclination and low black-hole mass) contribute. An improvement in 
the quality and quantity of the optical spectroscopic data could allow determi- 
nation of the H/9 transfer function, which could allow discrimination between 
these competing models. 

The differences between the characteristics of the H/? emission line on one 
hand and of the high-ionization lines such as He II A4686 on the other suggests 
a two-component BLR, which has been proposed on numerous occasions on 
other grounds (|2|). In this particular case, an interpretation that is at least 
qualitatively consistent with all the data and relatively simple is that the 
Balmer lines arise primarily in material that is in a flattened disk-like config- 
uration at a low inclination (to account for the narrow width of the H/5 line) , 
and the high-ionization lines arise in an outflowing polar wind, of which we see 
preferentially the component on the near side of the disk (to account for high 
velocity and blueward asymmetry) (^]). If indeed the Hf3 emission arises pri- 
marily in a low-inclination disk, our virial mass estimate of M = 1.1 x 10 6 M 
might seriously underestimate the black-hole mass. It could then be inferred 
that the NLS1 class might be best explained as low- inclination rather than 
low-mass, high accretion-rate systems. The strong rapid X-ray variability of 
NLSls seem to favor the latter explanation, but we note that at least some low- 
inclination accretion-disk models predict relatively strong, variable EUV/soft 



X-ray fluxes ([TT|; |6]), consistent with observations of NLSls. 



4 Summary 



On the basis of three years of combined X-ray and optical spectroscopic mon- 
itoring of the narrow-line Seyfert 1 galaxy NGC 4051, we reach the following 
conclusions: 

(1) The rapid and strong X-ray variations that characterize narrow-line Sey- 
fert 1 galaxies are detected in our X-ray observations of NGC 4051, but 
are not detected in the optical, consistent with previous findings (§). 

(2) On time scales of many weeks and longer, there does appear to be a 
correlation between the X-ray and optical continuum fluxes. 

(3) The variable part of the B./3 emission line has a Doppler width of about 
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1100 km s -1 , and a time delay relative to the continuum of about six 
days. Combining these quantities leads to a virial mass estimate of ap- 
proximately 1.1 x 10 6 M & . 

(4) The HellA4686 emission line is strongly variable, although an accurate 
time delay cannot be measured. This line is about five times as broad as 
the H.f3 line, and is strongly blueward asymmetric, as are the UV high- 
ionization lines in this object. 

(5) In the BLR radius-luminosity relationship, we find that narrow-line ob- 
jects (those with V FWUM £ 2000 km s -1 ) seem to fall on the same locus as 
AGNs with broad lines. 

(6) In the virial mass-luminosity relationship, narrow-line objects populate 
the low-mass end of a rather broad envelope; they have virial masses typi- 
cally an order of magnitude lower than other AGNs of similar luminosity. 

(7) During the third year of this program, both the X-rays and He II A4686 
nearly disappeared, while the optical continuum and broad H/5 emis- 
sion line were only slightly fainter than previously and continued to vary 
significantly. This suggests that the innermost part of the accretion disk 
went into an ADAF state, greatly reducing the production of high-energy 
continuum photons from the inner part of the accretion disk. 

A picture that is consistent with the emission-line characteristics is one in 
which the Balmer lines arise primarily in a disk-like configuration seen at 
low inclination (i.e., close to face-on), and the high- ionization lines arise pri- 
marily in an outflowing polar wind. The high-ionization lines are blueward 
asymmetric because we see emission preferentially from the near side, with 
the far side at least partially obscured by the disk component (which might 
be an extension of the accretion disk). If this interpretation is correct, the 
virial mass estimate of M — 1.1 x 1O 6 M is lower than the actual black hole 
mass by a factor of sini. If NGC 4051 is typical of the NLS1 class, then it 
might be that NLSls are best described as low-inclination rather than low- 
mass, high accretion-rate systems. While the observed X-ray characteristics of 
NLSls rather favor the latter explanation, they are not inconsistent with an 
inclination-based origin. Indeed, the full explanation of the phenomenon may 
involve both inclination and black-hole mass. 

For support of this work, we are grateful to the National Science Foundation 
(grant AST-9420080 to The Ohio State University). 
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